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Results are presented from a numerical simulation of transition control in plane 
channel and boundary layer flows. Details of the channel flow control are available 
in reference 1. The analysis is based on a pseudo-spectral/finite-difference semi- 
implicit solution procedure (ref. 2) employed to numerically integrate the time- 
dependent, three-dimensional, incompressible Navier-Stokes equations in a doubly 
periodic domain. In the channel flow, we find the active periodic suction/blowing 
method to be effective in controlling strongly three-dimensional disturbances. In 
the boundary layer, our preliminary analysis indicated that in the early stages, 
passive control by suction is as effective as active control to suppress instabili- 
ties. Our current work is focused on a detailed comparison of active and passive 
control by suction/blowing in the boundary layer. 
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GOVERNING EQUATIONS 


CONTINUITY EQUATION 

3u j _ 
oX i 

MOMENTUM EQUATION 

!!i + !!i= . i iL_ + A i 

at u j ax. p ax. v ax. ax. 

J ' J J 

THERMAL-ENERGY EQUATION 

3T il_ = JL_ 3 2 T 
at u i ax. pC ax. ax. 

' r J J 

CONSTANT PROPERTIES; NO VISCOSITY; TEMPERATURE FIELD UNCOUPLED 
EQUATIONS NONDIMENSIONALIZED BY U Q , h 

FLOW DRIVEN BY A CONSTANT MEAN PRESSURE GRADIENT 2/Re, Re = U Q h/v 
CONVECTIVE TERMS PUT INTO A FORM THAT CONSERVES ENERGY AND MOMENTUM 

9u i ^ , 9u i 9u j, . 3P* _2. I 1 9 ^ u i 

at u j 'ax. ax.' “ ax. Re il Re ax. ax. 

J I I J J 

P* = P/p + \ U jUj 
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SOLUTION PROCEDURE 


SAME TECHNIQUE AS THE VELOCITY FIELD 

Non-dimensional ized with (T-T )/(T w -T q ) ; h and U Q 
Adams-Bashforth 2-step method for the advective terms 

Crank-Nicholson implicit scheme on the diffusive terms. j 

Periodicity along allows x-j and j 

• Two-D Fourier transform in the x^-x^ plane j 

• The Pseudo spectral method in the x^j x^ directions 

• Finite differences with variable mesh along the direction ; 

Solution in Fourier space as a tridiagonal system I 

Back transformed into physical space to obtain temperature = 

field at (n+1 ) ] 

IMPLEMENTATION OF BOUNDARY CONDITIONS | 

No-sl ip B.C. 

U 1 = u 2 = u 3 = 0 
p from Xp-momentum equation 
. Suction B.C. 

Flow homogeneous along x^ , x^ 

• Incoming mass flow rate must equal to the outgoing mass flow rate 

3<u 9 > 

L_ = o 

ax 2 

<U 2 > = CONST, or g(Xj, x^) 

, Hence velocity magnitude and direction at one wall must be preserved 
throughout the flow field to satisfy continuity 

. Physically plausible condition is suction-blowing or periodic b. cond. 
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TRANSITION IN WALL-BOUNDED FLOWS 


2-D Tol lmien-Schl ichting waves. 

Formation of streamwise vortices. 

Formation of shear layers away from the wall due to 
vortici ty-induced velocity. 

Secondary instability 
(kinks and spikes). 

Breakdown into smaller scales, formation of wall shear, 
hairpin eddies. 

Turbulent spot - horseshoe vortex-turbulence. 


MODEL PROBLEMS 


* Periodic plane channel flow 


* 


Periodic boundary layer. 






MAXIMUM PLANE-AVERAGED RMS VELOCITIES 

The temporal development of plane- averaged maximum velocities is presented. 
These velocities provide comparisons between the controlled flow (three-dimensional 
control) and no- control cases for the channel flow. 



No control 

a. u 1 RMS 

b. u 3 RMS 

c. u 2 RMS 


With control 

d. u-j RMS 

e. u 3 RMS 

f. u 2 RMS 


c 
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TEMPORAL DEVELOPMENT OF U- COMPONENT FLUCTUATING VELOCITY AND ITS HARMONICS 

The evolution of the various Fourier modes indicated that all amplitudes are 
significantly reduced, and after the third control wave they all decay rapidly. 


No control 


With control 







NORMALIZED ONE-DIMENSIONALIZED WAVE SPECTRUM 


In the uncontrolled case we observe that energy transfer to the high wave num- 
bers is indicated by a full spectrum. In the controlled case, this does not occur 
and energy is concentrated in the low wave numbers preventing the development of 
higher harmonics. 


SPECTRUM OF u! VELOCITY * -0.9 



+ - ♦(K 1 )/+(l) 

a. Without control b. With control 
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DISTRIBUTION) 


plitting develops. The control wave 
reduces the amplitudes. The uncontrolled 


e naxinun ur rhs along X3 at t=50. 



X3 


O MAX IMUfl Ul* RMS ALONG X3 AT T=0O . 
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^-FLUCTUATIONS ALONG X ± 

In the uncontrolled case a strong negative spike develops between T = 40 and 
T = 50, No evidence of spike formation and nonlinear distortions is observed in the 
uncontrolled case. As T = 60, the controlled distribution is nearly sinusoidal, 
whereas the uncontrolled case shows a broad frequency content. 


FLUCTUATIONS ALONG XI AT T*20. 


U1 FLUCTUATIONS ALONG XI AT T = 5C 


FLUCTUATION ALONG XI T=40 


FLUCTUATIONS ALONG XI AT T=Q0 


With control 


With control 











SURFACES OF CONSTANT TEMPERATURE 


In these figures, three-dimensional representations of the temperature field 
(treated as a passive scalar) are displayed. The uncontrolled flow displays evidence 
of strong mixing and a highly convoluted temperature surface, while the controlled 
flow is relatively uniform and indicates local laminarization. 



( ) 4 7T (-)4t> 
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SURFACES OF CONSTANT TEMPERATURE (CONC . ) 


No control 


T = 40 


With control 







CONCLUDING REMARKS 


* 2D and 3D wave interaction in channel flow. 

* Preliminary calculations indicate comparable effects of passive and active 
control in the periodic boundary layer. 

* Passive temperature effective to tag flow dynamics. 

* Need for space-evolving numerical experiments. j 

! 
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